Absence of spontaneous magnetic order of lattice spins coupled to itinerant 
interacting electrons in one and two dimensions 
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We extend the Mermin- Wagner theorem to a system of lattice spins which are spin-coupled to 
itinerant and interacting charge carriers. We use the Bogoliubov inequality to rigorously prove that 
neither (anti-) ferromagnetic nor helical long-range order is possible in one and two dimensions at 
any finite temperature. Our proof applies to a wide class of models including any form of electron- 
electron and single-electron interactions that are independent of spin. In the presence of Rashba or 
Dresselhaus spin-orbit interactions (SOI) magnetic order is not excluded and intimately connected 
to equilibrium spin currents. However, in the special case when Rashba and Dresselhaus SOIs are 
tuned to be equal, magnetic order is excluded again. This opens up a new possibility to control 
magnetism electrically. 
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Since the seminal work on phase transitions by Ho- 
henberg [1] and Mermin and Wagner [2] it has be- 
come common knowledge that spontaneous order in low- 
dimensional systems is generically not possible at any 
finite temperature. In these studies, the use of the Bo- 
goliubov inequality [3] was essential: Hohenberg used it 
to rule out superfluidity [l]and Mermin and Wagner to 
rule out magnetic order in Heisenberg spin systems [2] in 
dimensions d < 3. This approach is very powerful and 
was then applied to many different systems [4-10], in- 
cluding the Anderson and Kondo lattice models [11, 12]. 

For systems in the continuum, the weak coupling ap- 
proximation is often applied leading to an effective ex- 
change coupling between the localized spins which is of 
the RKKY-type [13]. RKKY interactions occur in many 
physical systems, prominent examples of present interest 
are heavy- fcrmion systems [14], diluted magnetic semi- 
conductors [15-18], and nuclear spins in low-dimensional 
conducting nanostructures [19-21]. The latter system 
plays an important role as noise source for spin qubits 
in GaAs or InAs quantum dots [22-24], and much effort 
goes into understanding and controling the nuclear spin 
bath, with one possibility being to freeze out the nuclear 
noise by magnetic order [25, 26]. 

In contrast to the Heisenberg exchange, however, the 
RKKY interaction is long-ranged and thus is not covered 
by the original Mermin- Wagner theorem which requires 
the spin interactions to decay sufficiently fast with dis- 
tance r (faster than l/r 2+d ) [2]. Addressing precisely this 
issue, Bruno [10] was able to rule out in RKKY systems 
magnetic order in one dimension. A similar conclusion, 
however, for the two-dimensional counterpart appears 
still to be missing. Here we will fill this gap by rigorously 
proving the absence of order for a rather general class of 
systems which consist of lattice spins embedded in a con- 
tinuum of itinerant electrons with which they interact by 
an isotropic on-site spin interaction. The allowed electron 



Hamiltonian H e is very general and may include electron- 
electron interactions as well as any single-particle poten- 
tial (such as lattice or disorder potential) that does not 
depend on spin. For this class of models we prove then 
that in the thermodynamic limit fcrro- and antifcrromag- 
nctic, as well as helical, long-range order of the lattice 
spins is excluded at any finite temperature in dimensions 
one and two. We show that this conclusion remains valid 
when short-range Heisenberg interaction between lattice 
spins is included. Our result also applies to the RKKY 
case, since this regime is obtained from the full one by 
lowest order perturbation expansion in the on-site spin 
interaction [13] including the full H e [26]. 

Moreover, we consider the effect of Rashba [27] and 
Dresselhaus [28] spin-orbit interactions (SOI) which ex- 
plicitly break the spin symmetry. Our argument be- 
comes then inconclusive and magnetic order cannot be 
excluded. While this finding is not unexpected it is re- 
markable that it is closely linked to the existence of equi- 
librium spin currents studied recently in spintronics [29- 
31]. Even more remarkably, we find that in the special 
case when Rashba (a) and Dresselhaus (/?) SOIs become 
equal, magnetic order is excluded again. Since a can be 
electrically tuned to j3 [32-34] , this opens up a new way 
to tune magnetism by electrical gates. 

Finally, we note that the absence of spontaneous order 
proven here is valid only in the thermodynamic limit; 
thus, effective ordering in nanostructures of finite size at 
sufficiently low (but finite) temperatures is not in conflict 
with our findings. 

Model. We consider a lattice {Rj}^ filled with JVj 
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H) located at the sites Rj. The lat- 
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tice is embedded into a volume £1 containing N e itiner- 
ant electrons which couple to the lattice spins via on-site 
spin-spin interactions. The Hamiltonian for the entire 



system 


i reads, 












H = 


H e 


N, 


■ ij + h 




-iQR; 


7/ + h.c. 


), 


where H e 


= H Q + V + U = 




pf/2m + Ef<. 


,n 



(1) 



X^Ji C7(f i) is the Hamiltonian describing the electron 
system. Here, m is the mass and pi the momentum op- 
erator of the i th electron, Vij = V(ii — Yj) the electron- 
electron interaction of electrons at positions fj and Tj, 
and U(r t ) an arbitrary spin-independent single-electron 
potential. Typical examples for C/(f"j) are periodic lat- 
tice potentials, disorder potentials, electron-phonon in- 
teractions [35], etc. We remark that in contrast to pre- 
vious work on lattice models [11, 12], wc do not re- 
strict the motion of the electrons to the sites of a lat- 
tice (tight binding limit) but allow them to move in the 
real space continuum. Further, J denotes the coupling 
strength of the isotropic spin interaction at lattice site 
Rj, Hj = JEfiiSj • Ij, where Sj = S(Rj) is the 

electron spin density operator S(r) = X^i=i ^i<5(r — r,), 
with s, = (sf,sf,sf) being the spin-1/2 of the i th elec- 
tron. The vector components of each spin, s k and /j-, 
satisfy standard spin commutation relations. Finally, 
to probe the order for the lattice spins Ij we break 
the symmetry by an external (fictitious) field h point- 
ing in, say, the z direction, which we let then go to 
zero at the end. This leads to an additional Zeeman 
term H Z (Q) = hj^li e'^^I* + h.c. To rule out fer- 
romagnetic order we will choose Q = 0, whereas to ex- 
clude antiferromagnetic order we will choose Q such that 



-iQR 



-1, if R connects sites from the same sublat- 



tice, and e ~ 4< 3 R = — 1, if R connects sites from different 
sublattices. 

To prove the absence of spontaneous order for the lat- 
tice spins Ij we follow Ref. [2] and make use of the Bo- 
goliubov inequality [3] , which is an exact relation between 
two operators A, C, and a Hamiltonian H, 

\{{A,rf}){[[C,H]^]) >k B T\{[C,A))\\ (2) 

Here, (A) = Tre- H / ksT A/Tre- H / ksT denotes the ex- 
pectation value in a canonical ensemble, T the tem- 
perature, Ub the Boltzmann constant, and {A, B} = 
AB + BA the anticommutator and [A, B] = AB-BA the 
commutator. It is assumed that all expectation values are 
well-defined and exist in the thermodynamic limit defined 
by N e , Nj, Q — > oo with electron density n e — N e /tt and 
density of lattice spins rii = Ni /CI finite. 

Proof- The strategy of the proof consists of using the 
Bogoliubov inequality to derive an upper bound for the 
order parameter corresponding to the phase transition we 
want to discuss. If this bound turns out to be in contra- 
diction with the presence of long-range magnetic order, 
then the absence of the corresponding phase transition is 



rigorously demonstrated. The success of the procedure 
depends crucially on the choice of the operators A and 
C in (2). As we shall see, the appropriate choice for our 
case is given by 

G q = D_ q + i_ q + D_ q + i_ q , Aq = J q+ Q + -'q-Qj (3) 

where the Fourier transforms are given by S q — 
£^i e -iq-fig. and i q = j-J2 i e-ii-njf. [gg]^ and where 

B ± = B x ±iB v . Note that C q and Aq are not Hcrmitian 
in general. Since the Bogoliubov inequality (2) is valid 
for any wave vector q, it can be generalized to 



^({A q ,4})>fc B T 
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<[[C q ,tf],C q ] 



(4) 



where the sum runs over all q's in the first Brillouin 
zone of the reciprocal lattice. We note that the above 
choice for C q and Aq is essential also for the following 
reason. Besides the fact that ^2„([Cq, A q ]) can be ex- 
pressed in terms of the lattice spin magnetization, the 
generally complicated interaction terms V and U in H e 
simply drop out of the calculation since they commute 



with C, 



q> 



[C q , H e ] — [5_ q + S_ q , Hq\. 



(5) 



This simplification is a crucial advantage of first over sec- 
ond quantization formalism since spin and position oper- 
ators of the electrons trivially commute. (Note, however, 
that the expectation values still contain the full Hamilto- 
nian including U and V.) Hence, our proof goes through 
for any form of the potentials V and U as long as they 
are spin independent. 

We now focus on the various terms in Eq. (4) and 
find bounds for them. Here, we outline only the main 
steps of the calculations and defer details to the Ap- 
pendix [37]. As a first step, let us evaluate the dou- 
ble commutator on the right-hand side of inequality 
(4). By virtue of the commutation relation [S_„,Hq] = 
-2^E 4 ^{p*,e"i-H we obtain that [[C q , H e ], C q ] = 

— N e q 2 . The part of the double commutator with Hj 
vanishes since [Cq,Hj] = 0. Indeed, [S_ q ,iJj] — 
»Ei J -e <q -* i 5(ri - Rj)(Ij x s 4 ) ± , and thus [S± q ,Fj] = 

— [I_ q ,.Hy]. After some calculations (sec [37]) wc find 



that [[C q ,F,(Q)],C q ] = 
Hence, 

<[[C q ,ff],Ct]) = JV e 



-4M£j e 



-iQRjjz + hc _y 



I g 2_ 4/l ^ wf(Q) 



(6) 



where the lattice spin magnetization appearing in Eq. (6) , 
which wc identify as the order parameter, is defined by 
"if (Q) = wj(T,j e-^^i* + e iQ ' R ' /?). The commuta- 
tor on the right-hand side of inequality (4) can also be 
expressed in terms of mf (Q), 



([C q ,i4 q ]) = -2JV J mf(Q). 



(7) 



Finally, the sum on the left-hand side of Eq. (4) can be 
bounded as follows, 

^({^ q ,4}) = 2iV / ^({7+,J-}(l + co S (Q-R j ))) 



< AN]{2lf 



(8) 



where we have used that £ e* q '^ R * Rj ' = A^r^r., 

and ({If ,17}) < (2I) 2 . Using Eqs. (6), (7), and (8), we 
obtain from the Bogoliubov inequality (4) 

w?(2i) 2 /2 > k B ry 4J V/™/(Q) 2 , (9) 

Our goal is to rule out spontaneous magnetization in the 
lattice spin system, therefore we are interested in the 
behavior of the order parameter mf (Q) in the limit of 
vanishing external field, i.e., h — ► 0, after we have taken 
the thermodynamic limit. We need to distinguish two 
cases: i) mf(Q) = 0, V/i around h = 0; m) mf(Q) 7^ 0, 
V/i around /i = 0. If i) is satisfied, there is no order 
and the proof is completed. If ii) is satisfied, we need to 
show that lim/ l _ 5 .o mf (Q) = follows from inequality (9) 
in the thermodynamic limit. In this limit, the sum can 
be replaced by an integral, 



[2lf > 



kBTN lV 
N e (2w) 



d d q 



mf(Q) 2 



|q|<|qd 



2m 



|zmmf(Q)| 



(10) 



where v = 2Nj/N e , q c is an arbitrary cut-off vector lying 
in the first Brillouin zone, v — Q,/Nj, and we have used 
that ([[C q ,i7],C q ]) < N e (q 2 /m+ \2vhm\ (Q)|). In the 
one-dimensional case (d= 1), Eq. (10) gives 



\isj\h\ 
T 



arctan 



v/2H 



> 



im? 



^Pi' 



|Wimf(Q)| 

where Ai = iv(2I) 2 n e y/v / '(ksV^rn) . In the limit h — >• 0, 
the left-hand side of inequality (11) vanishes and this 
implies that lim/j_>.o f/ (Q) = 0. The two-dimensional 
case can be treated in a similar way. For d = 2, inequality 
(10) leads to the following relation 



A2 
T 



log 



1 



2m|i / /imf (Q)| 



>m?(Q) 2 , (12) 



where A2 = 2\f2\\/ ^Jvm. It follows from inequality (12) 
that lini/j^o m f(Q) = here, too. Since our arguments 
were independent of the choice of Q, we have proven 
that neither ferromagnetic nor antifcrromagnetic long- 
range order of the lattice spins is possible at any finite 
temperature T > in one and two dimensions. 

The absence of order can be traced back to the in- 
creased fluctuations in the lattice spin system in lower 
dimensions. These fluctuations, in turn, have their origin 
in the kinetic energy of the electrons, as one can explicitly 



see from Eq. (10) where the term q 2 /2m is responsible 
for the divergency in above q-integrals for d = 1 and 2. 

Next, we show that helical long-range order of the 
lattice spins is also excluded. The strategy of the 
proof remains the same and we shall be brief (for de- 
tails see [37]). To study this type of order, we con- 
sider the symmetry breaking Zeeman term Hz(Q) = 
y/2/3h^2 ■ e~ l( 3' Rj / + + h.c. and the magnetic order pa- 



iQRj f+ 

j 



which 



rameter mf(Q) = 

corresponds to a spin helix in the a;y-plane. Note that the 
spin part of Hamiltonian (1) is isotropic and consequently 
all choices for the helix arc equivalent. The operators C q 
and Aq for the Bogoliubov inequality (4) are now chosen 
to be 

■q and ^q = -7- fc+Q - Vq) • ( 13 ) 



Cq = Sl 



It 



The double commutator on the right-hand 
side of Eq. (4) becomes then ([[C q ,#], C q ]) = 

N e (q 2 /4m-vhm I (Q)/2). Since ([<5 q ,! q ]) 

(Nr/V2)mj(Q) and £ q ({l q ,I q _}) _ < 2N 2 (2I) 2 , 
Eq. (4) takes in the thermodynamic limit exactly the 
same form as Eq. (10), where mf(Q) must be replaced 
by mj"(Q). We thus conclude that lim/,^0 rnj~ (Q) = 
for any Q and hence long-range helical order is also 
excluded in one and two dimensions at any T > [38]. 

As a further generalization, short-range impurity-spin 
Heisenberg interaction Hz = Y]j jlj-ji-j ■ Ij is added 
to Hamiltonian (1). When the couplings 2^ satisfy 
1/iVj £ t - ■ \Xij I (Rj — Rj) 2 < 00, then both proofs to 
exclude (anti-) ferromagnetic and helical ordering re- 
main valid and lead to Eq. (10) with renormalized mass 
m* = ro/(l + 8ml 2 ^^ £„ %j\{^i - ^jY) [37]. 

Presence of spin-orbit interaction. Next we investi- 
gate the question of magnetic order in a low-dimensional 
electron gas in the presence of Rashba [27] and/or Dres- 
sclhaus [28] spin-orbit interaction which break the ro- 
tational spin symmetry of the Hamiltonian (1) explic- 
itly. The spin-orbit Hamiltonian is given by -ffso = 
H R + Hj}, with H R = aE?=i(p y i$i ~ Pfs v i), H B = 



•,N t! 



p\s\), where a (/?) is the Rashba (Drcs- 
for C q , we ob- 



0£r=i(pfs? 

sclhaus) coefficient. Using Eq. (3) 
tain [[C q ,i?so],C q ] = 4ma]J =0]I + 4m^ =0y , where 
we have defined the spin-current density operator as 
j"( r ) = ~\i Ya=\ sf{pi, 5(ri - r)} and its corresponding 
Fourier component j q = ^ £, sf{Pi, e _lqri }. These 
spin currents may lead to an intrinsic cut-off for the fluc- 
tuations in q, and thus help to establish order. To see 
this, we evaluate now the spin currents perturbatively 
around the free electron limit, i.e. U,V, J — 0, and at 
T = [391, 



1 -v \ r, mE F 
Oq=o,x)o = V— — a 



{jI= 



1=0, y I 



4-7T 
47T 



(14) 
(15) 



where Ep is the Fermi energy and the results are valid 
in the regime ma ,m(3 2 <C h 2 Ep [40]. Performing now 
a perturbative expansion in the parameters V, U, J, T 
around above free case, we conclude that (J*) ^ and 
(Jx) 7^ [41]. (In passing we note that in the station- 
ary and homogeneous limit, the spin-currents satisfy the 
relations (j%) = — (jy) and (j%) = — (jy) due to a general- 
ized continuity equation; see [37].) As a consequence, the 
commutator ([[C q , i?so],Cq]) appearing in Eq. (9) docs 
not vanish anymore and thus provides an intrinsic cut-off 
to the q-integral [cf. Equation (10)]. Hence, the bound 
for the order parameter we extract from inequality (10) 
is a constant which does not vanish in the limit h — > 0. 
Thus, our argument becomes inconclusive and we cannot 
rule out (anti-) ferromagnetic order in this case. 

Similarly, for helical order our argument remains in- 
conclusive, since [[Cq,H so ],Cl] = ma(j v ti=0tX -j^ =0ty ) + 

m l3(jq=o y ~ 3q=o.x)i which, will not vanish in general. 

Next, let us consider the special case a = (3 where 
new symmetries emerge [42]. Then, the leading terms, 
Eqs. (14), (15), cancel, indicating that the physics 
changes dramatically. Indeed, by making use of the 
"gauge transformation" U = e 4 ^*= Afc ' rfc , where A& — 
— am(s% — §f)(l, 1, 0), to remove the SOI from the Hamil- 
tonian, we can prove as before [37] that (anti-) ferro- 
magnetic order in the z direction can now be excluded 
rigorously for any T > and d — 1,2. Similarly, 
we can rule out helical ordering described by the or- 
der parameter mf = ■^j-(^2 1 e~ lCi ' Ri lt + h.c.) with 
Q = ^/2am(l, 1,0) [for rotated coordinates (x,y,z) —¥ 
(x', y', z>) = (z, (x + y)/V2, (x - y)/y/2); see [37]]. 

Thus, quite remarkably, this spin-orbit effect suggests 
the control of magnetism by electrical gates, namely, by 
tuning the Rashba SOI (a) [32-34] from the regime a^/5 
(ordering not excluded) to a — (3 (ordering excluded). 

Conclusions. We proved an extension of the Mcrmin- 
Wagner theorem for lattice spins interacting with itiner- 
ant electrons, and showed that spontaneous order of the 
lattice spins is ruled out in one and two dimensions at 
finite temperature. In the presence of Rashba (a) and 
Dresselhaus (f3) spin-orbit interactions, however, sponta- 
neous order could not be excluded, unless for a = f3. 
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I. MODEL HAMILTONIAN 

To make the Appendix largely self-contained we restate the problem defined in the main text briefly. We consider a 
lattice {Rj}j=i filled with Ni spins L, = (H, I? , H) located at the sites Rj. The lattice is embedded into a volume fi 
containing N e itinerant electrons which couple to the lattice spins via on-site spin-spin interactions. The Hamiltonian 
for the entire system reads, 

N, 

H = H e + H I + jJ2Srh+ H z(Q), (1) 

i=i 

where H e = H a + V + U — J2i=i P?/2w + J2i<j V%j + J2i U{ii) is the interacting electron gas Hamiltonian, Vij = 
V(ii — ij) the potential describing the interaction between two electrons at position f j and Tj, respectively, and U(ii) 
is an arbitrary potential for an electron at position rj. Further, J denotes the coupling strength of the isotropic 
spin interaction at lattice site Rj, Hj = JYli=i Sj ■ Ij, where Sj = S(Rj) is the electron spin density operator 
S( r ) = Sj=i ^i^( r — *i) with Sj = (sf,sf,sf) being the spin-1/2 of the i th electron. The lattice spin Hamiltonian 
Hz = Yljj 1-ifi-j • Ij describes an isotropic interaction between lattice spins with coupling constants Z^ satisfying 
1/Nj J^. . |Ijj|(Rj — Rj) 2 < oo. Finally, the Zeeman term Hz(Q), accounts for the presence of an external (fictitious) 
magnetic field which breaks the symmetry of the lattice spins. To rule out (anti-)ferromagnetic order we choose 
H Z (Q) = ft£j2i e- iQR *I? + h.c. and to cxlcude helical order we choose H Z (Q) = h£j2i e~ iCi - n df + h.c. 

II. BOGOLIUBOV INEQUALITY 

For the proof to rule out order we follow Ref. [1] and make use of the Bogoliubov inequality [2], which is an exact 
relation between two operators A, C ', and a Hamiltonian H, 

l -{{A,ti}){[[C,H],C>]) >k B T\{[C,A])\\ (2) 

where, (A) = Tre~ H/>kBT ' A/Tre~ H / knT denotes the expectation value in a canonical ensemble, T the temperature, ks 
the Boltzmann constant, and {A, B} = AB + BA the anticommutator and [A, B] — AB — BA the commutator. It is 
assumed that all expectation values are well-defined and exist in the thermodynamic limit defined by N e , Ni, fl — > oo 
with electron density n e = N e /Vt and density of lattice spins nj = iVj/fi finite. 

If the operators A and B depend on the reciprocal vector q, then the Bogoliubov inequality can be generalized to 

where the sum runs over all q's in the first Brillouin zone of the reciprocal lattice. 

III. (ANTI-)FERROMAGNETIC ORDERING 

In this case, the Zeeman term is given by Hz(Q) = h^j e~ 4( 3' Rj ' I? + h.c. and the magnetization is defined by 

™ f (Q) = W \ ^ ( e " Q ' M ! + e * Q ' R ^) ) • (4) 



We choose C q and Aq in the Bogoliubov inequality (3) such that 

Cq = 5_ q + i_ q + S_ q + I_ q and A q = I q+ q + ^ q -Q- (5) 

Oq Oq 

The strategy of the proof consists in evaluating every term which enters the Bogoliubov inequality (3) to derive an 
upper bound for the order parameter corresponding to the phase transition we want to discuss. If this bound turns 
out to be in contradiction with the presence of long-range magnetic ordering, then the absence of the corresponding 
magnetic phase transition is rigorously demonstrated. 

A. Evaluation of [[C* q , H], (Cg) 1 "] 

We decompose the double commutator [[C q ,-ff], (C q )^] into the following four parts which we calculate below, 
i)[[C-+C+,H e ],(C-+C+)*], ii)[[C-+C+,H A ],(C-+C+?], m)[[C-+C+,ffz(Q)],(C-+C+)t], and (6) 

iv)[[C-+C+,H x ]. (7) 

1. Evaluation of i) 

Let us first determine [C~ , H e ] , 
[C~,H e ] = [Sl q + JZ q) H + V + U\ = [J2 e**$r + J] e^M", # + V + V\ = $2 e^s~, £ p 2 fe /2m] 

i j i k 

= E t^>^] = E£p- [e^.W +E £[^.W ' ft" ( 8 ) 

i i i 

Since 

[e lq - f %Pi]?A = (e iq -^ Pi - p^e 4 ^)^ = -ie l * u ^ - qe'^V + ie^ip' = -qe^V, (9) 

we can conclude that 

[C-,H e ] = -dY,^{Pi,e i ^}, (10) 

H ^-^ 2m 

and thus that 



£j I I I — III 

i i k 

= -J-T, ( ^si\q-{p i y^},^^} + [sr,s+}e-^ q -{p i y^} ) . (11) 

4 V A B I 

We can now evaluate expressions A and B, 

A = (l/2-«?)q- ([fce^'.e-****] + [e"^p t , e^*]) = (1/2 - J*)q • (e^'fre"***' - e^'fre***') 

= (1/2 - s*)q (Pi - Q- Pi - q) = -2(1/2 - s*)q 2 , (12) 

B = -2*fq • (e^'pie***' + Pi ) = -2sf (2q • p 4 + g 2 ). (13) 

s • 

Pi+q 

With the help of the above results, we can finally conclude that 

[[C-,H e ], (C q -)t] = - ± £(1/2 - st)(-2q*) ± £(-2*?)(2q • p, + 9 2 ) = ^ £ Sffr + ^g 2 . (14) 



Let us now calculate [[C+ , H e ], (C+)^]. Since 

[C+,H e ] = Ef^-Ie^'.ftl + ^Ie^.Pil-ft) =-q-E^{ft,e*»*«}, (15) 



2m ' ! ' J 2m / ^-^ 2m 

we thus have that 



[C* H e ], (C+)t] = -J- 2 sfs- [q • {p„ e^}, e-^] - -L £[«+ $r] e -*»-*'q • {p 4 , e**} 

i % 

l - 2(1/2 + S?)(-2« 2 ) - -J- 2 2S?(2q ■ fr + q 2 ) = ^-N e q 2 - |j E ^ ( 16 ) 



2m ^— ^ ' 2m ■<— ^ * l 2m ' 2m 

i 

Let us now calculate the mixed term [[C+ _ff e ], (C~) 1 ']. From 



2m 
we have that 



[C+,H e ] = -q-E^iP-^'"}' ( 17 ) 



l - E § t$tl<l ■ {Pi, ^}, e-^} - A- ^l^^+le-^-q ■ {f>i, e**} = 0. (18) 



2m ^ * * L " Lr ' J ' J 2m 

An exactly analogous calculation shows that [[C~ # e ], (Cq")T] = 0. This finally allows us to conclude that 

[[C q ,H],Ct] = liV eg 2 . (19) 

2. Evaluation of ii) 
Let us calculate [[C~ + C+,Hj], (C~ + C+ )*]. As a first step we evaluate 

[C q + C+ , Hj] = J[C~ + C+, 2 Sj ■ tj]- (20) 

3 

From 

[5Z q , 2 S d ■ Ij] = 2 e***<$7> E S **(** - R ^] = E^s-^^ - R,)s, • I,] 

= E Ee^5(f i -R j )/ J a [5-,^]= E Ee^^-R^^^^e^-ie,^) 

a— x,y,z i,j a—x 7 y.z i,j 

= ^e'^^-R^iftxs,)-, (21) 

and 

[jr q ,EM;] - E eiq ' Rii r>E^(^- R ^)y= E E eiq ' Rj Wi-R;)[/7,/j] 

= E E e 4q ' R ^(^ - Rj)iS?i?(e xa( , - ie ya0 ) = E e iq - R ^(f 4 - R,)i(s, x I,-)", (22) 

a— x,y,z ij i,j 

we have that [<§!„, X) 7 - Sj ■ Ij] = — [-^Zq, X), Sj • Ij] and consequently [C~, Hj] = 0. A similar calculation for [C£,Hj] 
shows that 

[C+,Hj] = jEe iq ' fi *(ri-Rj)i(Ij xs,)+ + J^e'^^-R^jfsi xi 3 )+=0. (23) 

From this it follows that [[C~ + C+,Hj], (C~ + C+)t] = 0. 



3. Evaluation of Hi) 
Let us calculate [[C q , iJ z (Q)],C q ]. From 
[C q ,ffz(Q)] = [/I q + /± q ,#z(Q)] 



= hj^ (e iR J-( q - Q )(/- - /+) + e ^-(q+Q)(/- _ /+)) f ( 24 ) 



we obtain 

[[C q ,# z (Q)],C q ] = h[^ (>V(q-Q)(i- _ /+) + e ^-(q+Q)(/- _ /+)) ^ e -iqR, 7 + + g-iq-R*/- 

3 I 

j j i j 

= 2/t J] e -' QR ^ [/- , /+] + 2/i ^ e lQR < [IT , J+] = -Ah J2 (e _iQRj J/ + e iQR ' J?) . (25) 

J 3 3 

4- Evaluation of iv) 
Let us calculate [{Cz + C+,H X ], (C~ + C+) f ]. Since 

[C~ + C+,H T ] = [/+ q + t v J2 ZflI?I?] = [2 E ^'^ E ^'Wl 

= 2 2 ^e' q - R n4 x ,^n=2 XJ ^e jq - Rfc (K,W+^K,^]) 

a,k,j,l a,k,j,l 

= 2 X Ijie***' ie xa pifl? + 2 ]T ^e iq ' R < Ifie xa ^ 

a.j,l ce,j^l 

= 2 X ^;(e iq ' Ri - e iq ' R )(I; x !,-)*, (26) 



we obtain 



[[C q , Fx],C q ] - [2 X ^/(e iq ' Rj - e iq ' Ri )(t x *,•)*, 2]T e-*»- R * /£] 
= -4iX^e _iq ' R,c (e' qRi - e iq - R, )fe & x I,-),] 

k,j,l 

= -4iX^e- iqR He' qR ^ -e' q ' R ) fe 6 -^^' W + / "^' / /]> 
= -4i X I J -je- <q - R ' (e Jq ' R ^ - e iqR < )i £ ^pe^lf 

3,1 a 

-Ai J2 Ijie-**** (e iq - R ^ - e' qR * )* £ ^^e x ^I] 

3,1 a 

= -Ai X X,-ie- <q - R| (e* q ' R ^ - e jqR )i(/f /? + 7f/j / ) 



+4iX 2 "^ e_ * q ' Rj ( eiq ' Rj - e' q - R! )i(/f/| + Jf /?) 



4^^(e 4q,(R ^ R,) - l)(/f/| + Iffy) -4j]ljl(l - e iq -( R, - R ^)(Jf /? + If /J) 

-8]T^[1 - cos(q • (R, - R,))](7f J/ + /f/J). (27) 



The expectation value is then given by 

■,2,D . __ TJ. ,-> 
"2 



([[C q ,Fx],Ct]) = -8^X j7 (l-co S (q.(R 3 -R0))(/fi;+ifij)<8^|X j7 | g(Rj 2 R ° 2/ 2 

= 8/V £ 11,-jKRj - R) 2 = 8I 2 q 2 N! -±- ]T \l 3l \(K, - R) 2 . (28) 






<oo 



B. Evaluation of [C q , ^4 q ] 



Let us calculate [C q , -A q ]: 

[C v A,} = [/I q + /+ q , 7+ +Q + 7+_ Q ] = [7 Iq , 7+ +Q + 7+_ Q ] = £ (e-^ R ^ [7" , 7+] + e *Q** [7" , ^ 

= -2^ (e -iQ ' B ' J? + e iQ 'M?) . (29) 



C. Evaluation of E q {^q> ^U 
As a first step let us calculate {Aq, A„}, 



{A q ,A q } - I J2 (e- j(q+Q) ' R 'i+ + e -^-Q)- R '7+J ,£ ( e *(q+Q)-Mr +e *(q-Q)-M: 

= \^ /' e »(q+Q)-(R i -R,)/J+ ) J 7 \ + e «q-(R,-R ( ) e -iQ-(R+Ri){j : + ; f- j. X 



+ ^ ^ e iq-(R j -R ! ) e iQ-(R ! +R,){/+ ) /-} + e i(q-Q)(R,-R0{/+ ) 7-^ ( 30 ) 

With the use of £ qe *q-(R*-Ri) = N I 5 Rl ^ j , we obtain 

^{A q ,A q } = JV J [22{^,^} + E( e " i2Q ' Bi+ei2Q ' Rj K / i h ' i '7}) =^/E^' / 7K 2 + 2cos (Q- R J -))> 

q \ i j / j 

(31) 

and consequently 

/^{A q ,A q }\ < 4iV 7 £({/+, /-})<4iV 2 (2/) 2 , (32) 

\ q / i 

where we have used ({7>,7r}) < (2/) 2 . Indeed, {7r,7>} = 2(7f) 2 + 2(ff) 2 < 4/ 2 . 



D. Putting everything together 

The Bogoliubov inequality (3) takes now the following form, 

4iV 2 (2/) 2 /2 > k B TY^ = 4JV 7 2 mf(Q) 2 

In the thermodynamic limit we can replace the sum by an integral and obtain 

4^(2/) 2 /2 > fc f^ f dU N AN >^ 2 , (34) 

( 2 ^) d i|q|<|q c | ^<Z 2 -4WVimf(Q)' 

where v = fi/JVj, <i is the dimensionality of the system, m* = m/(l + 8I 2 Ni/N e m^- J^- |Zjj|(Rj — Rj) 2 ), and q c is an 
arbitrary cut-off vector lying in the first Brillouin zone. Since ([[C q ,_ff], C q ]) < N e (q 2 /m* + \2vhmj(Q)\), inequality 
(34) can be simplified to 



N e (2*) d 7| q |< M JtL + \vhmftQ 

where v= 2NijN e . 

IV. HELICAL ORDERING 

In this case the symmetry-breaking Zeeman term is given by Hz{Q) = V^AE, ( e~ lC *' Rj I^~ + e l( ^' Rj 1^ ) and 
the order parameter is a helix in the xy-plane: mf(Q) = -y/2/3-^- <5Z 7 - (e _lC *' R:, ij + e^'^/. - ) ). The choice of 
operators A q and C q in the Bogoliubov inequality (3) is as follows, 

f 



C q = 5i q + Ji q and A q = -^ (J+ +Q - 7 q _ Q J . (36) 



The numerical prcfactor l/v3 is chosen for later convenience. 

A. Evaluation of [[Cq, H e ], C?*] 
Let us first calculate 

[C q ,ff e ] = [^e^s^^p 2 /2m]=^|J-([e^,p,].p J +p l .[ e ^,p J ]). (37) 

i k i 

Since [e 1 *^ , p;] = -qe 4C i '*< , it follows that 

[C q ,ff e ] - -qE^{P,e^}, (38) 

i 

and consequently that 

[[5 q ,F e ],C q ] = -^[^s| q .{p i ,e^},^e-^lf] = -^^sf S f[q-{p i ,e i ^},e-^] 

i i i 

-is: ifis e "^ q - fe ' eiq '" } = -^? i ( - V) = i^ 2 - (39) 



B. Evaluation of [[C* q , flz(Q)], C q ] 



From 



we obtain 



[c q , ^z(q)] = E e iq ' Rj i;, v^73/iX) e ~' Q ' R ^/ + Vy3hj2 e * Q ' Rj ii} 

3 3 3 

= ^2/3h E e*< q - Q >- R ' If - ^/2jZh J2 e 4(q+Q) - R ^ Ij , 



[[C q ,-ff z (Q)],C£] 



/ V3 /i Zl eJ(q ~ Q) ' R "^ + - v^^E ei(q+Q) ' R 'f 'E e ~ iq ' R3 ^ 

i ! j 

V^/i E e -<Q ' R ' // - VW^ E e ' QR3 '7 • 



(40) 



(41) 



As a first step let us calculate 



C. Evaluation of [[C q , tfj], C* q ] 



[C q ,F J ] = [5i q + 7i q ,J^S j -I j ]. 



Since 



[s^^Sj-ij] = Q2e u **<§* i ,'£s k 5(h-& j )tj\= E Ee^ft-R,-)^,*?] 



3,k 



a=x,y,z i, j 



= E E e lq "^(r 4 - Rj)*#*f Cza/9 = E e***'^ - Rj)*(Ij X §,), 



a— x,y,z i,j 



and 



[/i^ESi-y = E eJq ' R ^E § ^- R ^)y= E E e<q ' R4 *('*- R j) i5 " f / c *^ 

i l id a—x 7 y,z i.j 

= E eiq ' R ^-R^(s 4 xi,) z , 
*,j 

we have [51 q , J2j Sj ' Ij] + [-^-qi Sj Sj • Ij] = and consequently [C q , Hj] = 0. 

D. Evaluation of [[C q , #z], <5 q ] 



From 



we obtain 



[C q ,fir] = E e i "- R *2 3 - J [^,^/n= E ^e iq ' Rfe ([4 2 ,/;]/f + /«[/|, Jf]' 

a,k,j,l a.k,j,l 

= Y^^^^pip? + J2 I » ei *' B4 *? ie *°e*? = *E W q ' Rj - eJq ' Ri )(ii x ij 

a,j,i <*,j,/ j,' 

[[C q ,Fx],C q ] = -,E^^ q ' Rfe (e Jq - R ^e 2q - R 0K,(i ; xi,) z ] 



fc,j,; 



-i E^ze-^ R ' (e* q ^ - e iq ' R ) E *W W^; 



-*E^ e ~' q ' Rj ( eJqRj - ^• R, )l]^e^ 1 /r/ 

jZ a 

-2E^(1 - cos(q • (R, - R)))(/JIf + /J/f )• 

3 1 



(42) 



(43) 



(44) 



(45) 



(46) 



The expectation value is then given by 

,2/"D . Ti.\2 



<[[C q ,ffz],C£]) = -2^Z j; (l-co S (q.(R J -R0))((ij : if + ijif))<2^|Z ji | g(Rj 2 R ' } 2/ 2 

= 2 ? 2 / 2 ]T |T j7 |(R,- - R,) 2 = 2 g 2 iY 7 / 2 -l ^ IZ^KB,. - R,) 2 . (47) 



E. Evaluation of [C q ,j4 q ] 
Let us now calculate [C q , Aq], 

[c q ,i q ] = Ee^7f,-^x) e "' (q+Q) ' R '^ - ^E e " i(, " q) " Rj ^i = 4E e ~ iQR ''/ + -^E elQRj v 

! 3 3 3 3 

(48) 
F. Evaluation of ^ {Aj,^,} 



Let us first calculate {yl q , A q }, 

^(/q+Q-Vqy ^ 



q'^q 



1 [ y^ e -i(q+Q)-Ri/+_y^ e -i(q-Q)-Rij- j ( J_ [ y^q+Cfl-Rj j- _ V^ e »(q-Q)-Rj f+ ] I 

= l y^ e '(q+Q)-(Rj-Ri){j+ ) /-}_ Iy^ e *q(R 3 -Ri) e -*Q(R J +R(){/+ /+} 
3 U 3 J,i 

_l^ e iq.(R J -R,) e iQ.(B,+R i )|/ j - j /-j + ly^ e i(q-Q).(R i -R,)|j' j - ) ^.j - (49) 

1,3 1,3 

By using ^ e *q-(Rj-Ri) = Ni6r 1} b.j, we obtain 

q \ 3 3 3 3 / 

\ J 3 3 / 

\ 3 3 3 J 

< ljv> (2E{/;,/-} + 4E(( J rj) 2 -( J r-) 2 ) + 4E{/;,/j'}) <2iv 2 (2/) 2 , (so) 



where we have used that 



{If, If} = 2{{I*f-{I]f± l {I*,I]}), (51) 

((I]) 2 - (If) 2 ) < 2/ 2 , (52) 

({/J, /J}> < 2/ 2 , (53) 

({//,/"}) < (2/) 2 . (54) 



G. Putting everything together 

With the use of the above results the Bogoliubov inequality (3) reads as follows, 



2N](2lf/2 > k B TY: ^Ti? 'laV (55) 



where to* is defined as above. In the thermodynamic limit we thus obtain 



k B TN lV f d mf(Q) 



2 



(2IY > , " / d d q— 2 ^^ , (56) 

where v and v are defined as above. 

V. PRESENCE OF SPIN-ORBIT INTERACTION 

In this section we want to investigate the effect of the presence of Rashba and Dresselhaus spin-orbit terms which 
break the rotational spin symmetry of Hamiltonian (1). The spin-orbit Hamiltonian under consideration is thus given 

by 

F SO = h r +h d = a j2(p^- m + p^mn - *ro> (57) 

i i 

where Hr and Hjy are the Rashba and the Dresselhaus spin-orbit Hamiltonian, respectively [3, 4]. 

A. ( Ant i-) ferromagnetic ordering 

In this case, the operator C q entering the Bogoliubov inequality (3) is given by C q = 5T q + IZ q + S-q + ^3 q - 
Let us first calculate [[C^Hr], (C q )^]. From 

[C q , H R ] = [5Z q + S± v H R ] = [£ e**{St + $,-)>«!>?«* - «)] 

i i 

= «E I (^ + KKle^M + [g + g,"^f] p?e^ - Qs+ + g t -)^[e iq "Spf] ~ [ff + ^ fflpf e'**< j , 

(58) 



and with the aid of 



(# + W - 2SfSf = i, (59) 

(# + $")*? = 2lf*J' = «? ) (60) 

[st + s-,%] = [25?,l?] = 2«f, (61) 



we have that 



[C q ,#*] = a^ (^[e il >'',p?]-«f[e^,^]-2tS?^e fa >-*') =«^ Q[e**#] -«?{p?,e*»*}) . (62) 
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Therefore, 



[[C q ,ff«],Ct] 



:Ef^ eiq '' < 'P?]-* 5 i^' eiq ' f< >)'E 



aE& T + s r)[^ ,q "%p?j,e 



o«q-ri ^S/l „-»qrj 



-Q 



E"*( S . t + 5 r)[{^.e <q '*'}.e- iq -*'] 



-2 9x 






2pf+fe 



<EW- 



(63) 



Let us now calculate [[C q , ifo], (C q )^]. Analogously to the above case we have that 

= ^E ( (^ + + *n^ q ' f %pf] +jg+ + gr,gf] ^e'^ - (5+ + s-)^[e iq - f %$] - [5+ + I", s\]p\e^ 



and consequently we obtain 

[[c^H D ],ci] = [/3E(^ eiqfs '^]"^ z {^' eiq '' s })'E e " q ' fs ( 

= ^E( 5 i' + *r)[^[e iq -* 4 ,^],e- i€, *«]-|9E"*( S i" + S . 7 )[{P?.c iq - fi },e- <q - f '] 



r -?+s: 



-2g» 






2p«+9 B 



(3J2^ 



We thus hnally conclude that 



[[C q) ffgo],Ct] = 4E(o^ + ^)S?. 



(64) 



(65) 



(66) 



B. Helical ordering 



In this case, C q is chosen to be C q = 5i q + 7f._. From 



[c q ,/f fl ] = [^^n^Y.inn -m)\ = «E a eiq ^>p^] - [e^i?,^?]) 



/ 



E 






[*?,*?] pfe** 



(67) 



we obtain 
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i ^ ' i 

El ■[■'•yf^.y iq-r, 1 — -iq-r,- -.si , -r A ^r,££ iqr, ~i — iqr, -.si 



a > -i 



\ A £33 



\ 



-2 9y 



2p;+q y 



+«E^I •^M{/';-''' qr 'K^' qr i + [«;'-^]'"' q - , M7>;^' qr } 



«E^-^i) 



-2g x 



2 Pf +9x 



Similarly, from 



[C q ,ff D ] = E^'^E (pfSf -OT)] = /j£ ([^'^Pf*?] - [e^f.M) 



/ 



stsne^ r \m + [strstm^ - m^^m - \sin}m* r 



= *E 



we obtain 



[[c q ,H D ],c q ] = [^E(^^' eiq " i } + ^?^' eiq ' f< })'E e_<q '' iS i] 

i ^ ' i 

= #E f^^{P^e^},e-^sf] + ^n$V q ' fi },e- iq ^| 



/ 



\ 



= />E 



a^s z r/«? *=«i- r 



V **? 



M,e'*'*},e- t * r <]+-[ 



P s z l p-^q-ri/^ P i q- r 



{jjf.e'*"} 



-2 9a; 



2pf+9x 



/ 






We can thus finally conclude that 



[[Cq,Hso],C J — —Hso- 



(68) 



(69) 



(70) 



(71) 



VI. CONTINUITY EQUATION FOR SPIN-CURRENTS 

Similar to Rcfs. [6, 7], we can derive a continuity equation for the total spin density operator S(r) = S(r) + I(r), 
where the lattice spin density is defined as I(r) = Si=i i?^( r ~ Rj)- Let us first recall that S a (r) satisfies the 
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Born-von Karman periodic boundary conditions and can thus be expanded as a Fourier series, 

^w = ^E e4q ' r ^> ( 72 ) 

q 

where S^^e-^sf. 

The Heisenberg equation of motion for the component a = x, y, z of the electron spin density is given by 

k a {v) = i[H,S a {v)]. (73) 

Since i[V,S a {r)\ = 0, i[U,S a {r)} = 0, and i[H z (Q), S a (r)] = 0, we only need to consider i[H ,S a (r)], i[Hj,S a (r)], 
and i[Hso,S a (r)}. Let us hrst calculate i[Ho,S a (r)]. By using the Fourier decomposition from Eq. (72) we obtain 

i[H ,S a (r)} = ^Y.^T.^l^ 1 ^- ( 74 ) 

q i 

From 

[p„ e -^]=-q e -^, (75) 

wc have 



i[H 0t S a {*)] = ^^^""ES^^'^' 6 ^" 1 



-4E e4q ' r( l-Jq = - V -J a ( r )' ( 76 ) 



ft 

q 



where the spin current density operator is defined as j Q (r) = X)i=i 2m^f iP*' ^(^* — r )} an d its Fourier decomposition 
reads j A (r) = jj X) q elq "Jqi wrtn Jq = 2mEt^{Pi' e ~ !qri } (^ = x iVi z )- We note that above definition of spin 
currents emerges naturally in the present context, and must be carefully distinguished from the one commonly used 
in the literature [5-7], j Q (r) = X^=i \{^h ^?^(^j ~ r )}' wnere the velocity operator, defined by Vj = i[H, r j]/h, picks 
up an additional anomalous spin-dependent term due to the presence of the SOI term i?so- This anomalous term is 
absent in j Q (r). 

Let us now calculate i[Hj, ^"(r)], 



(-7 

• £/3c«7 



i{Hj,S a (r)} = i[jY,hS(ik-n j )i j ,J2siS(h-r)} = J ^ E^ r * ~ ^j) 6 (Pi ~ r ) if j § l 

j.k i j3=x,y,z i,j 

= - j J2 E s ^ - r )^J e ^7 = jJ2 5 ^ - p )& A § j)«- ( 7? ) 

0=x,y,z j j 

Finally we consider i[Hso, S a (r)}. Since 

i[H R ,S a (r)} = ii^e^'Ifffl,^], (78) 

q 

and 

[H fl ,^] = aX)^^'E e " iq ' ffcS *]-°E^i'E e " <q '' fc5 *]' ( 79 ) 



with 



A = sfs«[pf , e^'* 4 ] + [If, *?]e-*»*# = ^ WZipf, e-<**<} (80) 

=*• «-e xa7 ^ §]{pl, e~ jqr *} = -ae xa7 2mj^ !/ , (81) 
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B = sVsflpt, e"***'] + [SI sf}e-^pf = ^ w7{pf, e"^} (82) 

=> \ ae y<*i E %&> e ~ 4q ' ?l } = ^<*evon2m% tX , (83) 

we have that 

i[H R ,S a (r)] = --ae xai 2mj^ + -ae yaj 2mj2, (84) 

where we have used j*(r) = ^ ^Z e * q r iq, 7 (7 = x ' 2/j z )- 

For the Dresselhaus term the calculation is similar. Since i[Hn, S a (r)] = i^ J2q eiq ' t [Hd, Sq] and 

[HD,§Z] = l352m%,e- i «- f *§?]-l3 , £\p i %,e- i «- f *S?], (85) 



=A =B 



with 



^4 = Sf S?K, e-***<] + [af, Sfle-^'p? = ^e xai s]{pf , e"**'*} (86) 

=> ^<Wy E «?{^. e _<q ' f4 } = ^W^,*' (87) 

B = SVjf^.e-^"**] + [^,^]e-^Sr = | W7{P?,e _iq *} (88) 

=> \fcv*i E *7{^r, e-^} = ^e BQ7 2mii tf , (89) 

i 

we can conclude that 

i[# D ,S Q (r)] = - l -Pe xai 2mjZ + l -(ie yai 2m 3 ^ (90) 

where we used j*(r) — jj ^Z e * qr iq, 7 as above. 

Let us now derive the Heisenberg equation of motion for the lattice spin density. Since i[H e , I a (r)] — 0, we only 
need to consider i[Hj,I a (r)] and i[Hz{Q,), I a (r)]. Let us first calculate i[Hj,I a (r)]: 

i[Hj,I a (r)} = t[jgM(fi-R J -)i J -,X)^( R i- r )] 

= iJ E E 5 ( ? « ~ R j)*( R i - r)isfi]e0 ai = -J E E 5 ( R J ~ v ) s p]^ a7 i3 

@=x,y,z i.j 0=x,y,z j 

= -jJ2S(R j -r)(i j xs j ) a = -i[Hj,S a (r)]. (91) 

For the (anti-)ferromagnetic case the Zeeman term is Hz(Q) = ft^j (e - ^ 1 ^/? + e l( ^' Rj 'J? J , and 

i[# z (Q), /«(r)] = »[h£(«" 9 'H+e <Q -H).E / W- r )] 

l 3 

= ihY,(e~ lQRl +e lQ ' R ')(5(R J -v)[If,If] 

3,1 

= »^(e- iQ ' R ' + e iQ - R 0*( R j - r)^a 7 /;. (92) 

The continuity equation for the z-componcnt of the total spin density takes the following form, 

E*(r) = -V • j z (r) - ae xzy mjy + ae yzx mj x - Pe xzy mj y x + Pe yzx mj y , (93) 

=> E*(r) + V • j*(r) = maUl + f x ) + raffl* + f y ). (94) 



14 

In the homogeneous and stationary limit, the left-hand side of Eq. (93) vanishes and this leads to 

<£> = -<#> and <#> = -<#>■ (95) 

For the helical case, the Zecman term is #z(Q) = \ZV3hJ2j (e~ iQ ' R: >lf + e iQR * 7j] and 

i[# z (Q),/ a (r)] - i[ v ^73/ l ^(e- i Q- R 'J+ + e^ R 'Jf),^/^(R j -r)] 

* i 

= zv^73^ ( e -^' R ^(R, -r)[i+,i«] +e iQ - R ^(R i -r)[Jr, j*]) 
j 

= Jv /273fe^ e- <Q - R ^(Rj - v){ie xai I] - e tfa7 /7) + ^v^/i^T e jQ ' R ^(R, - r)(ie xaj I] + e yaj IJ). 

3 J 

(96) 
For a = z we thus have 

i[#z(Q),/ a (r)] = z v ^73^( e -^ R ^(R,-r)(-/+) + e ^ R ^(R,-r)/-). (97) 



3 

In this case the continuity equation for the z-component of the total spin density reads 
±'(r) + V • j*(r) = ma(f y +f x ) + m/3(% + %) + iy/2/3hJ2 (e^^R; - r)(-/+) + e'^^R,- - r)/r) .(98) 

In the presence of the lattice spin-spin interaction term Hz, the right-hand side of Eqs. (93) and (98) acquires an 
additional term 

t [H x J a (r)} = i(£ J2 I#f,X;^(Rt-r)]=i £ VPRt-^tWl^l + I^W) 

j,l f3—x,y,z k j,l,k,[3 

= J2 2* W - r )^0jlfi? + J2 T H 5 i K o - y )^mIiI] = ^IjiSiRi - r)(Ij x /,)»■ (99) 

VII. EQUILIBRIUM SPIN-CURRENTS FOR [/, V, J, T = 

The aim of this section is to calculate uniform equilibrium spin currents for a system described by Hamiltonian 
(1) in the special case V, U,J,T = 0. The calculation follows the ones given in Refs. [5, 7]; however, due to the spin 
currents occcurring here without anomalous velocity term (see text after Eq. (76)), the leading term will turn out to 
be first order in the SOIs (instead of 3rd order [5, 7]). The single-particle Hamiltonian we consider takes the following 
form, 

H = h -r(a- x Py - VyPx) + T^yPy + <?xPx), (100) 

where p is the electron momentum operator and <J x ,y,z are the usual Pauli matrices. The eigenstates of the system 
are 

V M (r) - ^u s (k), (101) 

/ 1 \ /-f3k x + ak y \ 

where A is the area of the system and u s (k) = I r + ) > with r = -4= j3k y — ak x . The spectrum is given by 

\ s w^J v \ o J 

E s (k x , k y ) = %^+ sJ{a 2 + /3 2 )/c 2 - 4k x k y a/3, (102) 

2m, V 
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and it takes the following form in polar coordinates (k x ,k y ) = (kcos(9),ksm(9)) 



E s (k, 9) = —— + sk J (a 2 + B 2 ) - 4cos(6>) sm(9)aB . (103) 

2m v .. / 



-«(«) 



h 2 k±(ef 



The Fermi wavevectors k±{0) of the two branches are defined by Ep = E s=T i(k±(9)) = — ^-^ — h sk±(9)a(9) and 
the explicit expressions read 



, , „■. a[a)m / m ,„.\ Zm „ 

fe t W-±-y- + y(^aW) +^f- (104) 

In the perturbative limit where ma(9) 2 <C h 2 Ep, we expand k±(9) up to second order in a(0) and obtain 



**<« = V S^3^«< 8 » 2 * l"«" + V w E - (105) 

The spin-current tensor 7/ m is given by 

1 1 x— v /" <i 2 fc 

Tim = ^(jq=oj)o = ^ Z-, / 7^)2 ( a "» U * + ^ IJ m)k, SI (106) 

where (...)k,s denotes the expectation value in the cigenstates i/>k.s and the integration must be performed over 
k<k±(6). ' 

Let us first calculate T xy - Since v x — hk x /m, we have that {cr y ,v x } — 2hk x a y /m, and 

2ft , , . 2ft k T T v 2ft k T (Bk„ — ak T ) 

— ((Tykxh s = — s-krr = —s , ^- (107) 

With the use of polar coordinates we obtain 

271,,. 2ft /3/c 2 cos(<9) sin(6>) -afc 2 cos 2 (6») , . 

— ffi,% ks = — s . (108) 

m" m k V(a 2 + P 2 - 4cos(6») sin((9)a/3) 

S v • 

=a(0) 

From Eqs. (106,108) we obtain 

12ft 1 [ k - {6) [ 27r Bk cos(9)sm(9) - ak cos 2 (9) 



■"■" 2m(2^) 2 A +w y - a(0) 

ft 1 /" 27r ^/c_(6') 3 -A: + (0) 3 ^cos(6')sin(6')-acos 2 (l9) 
ac'- 



m (2tt) 2 7 3 a(0) 

ft 1 1 f 2n „ r 8m 3 ,„,, 12m 2 //1N1 /3cos(6»)sin(60-acos 2 ((9) , N 



m (2tt) 2 3 7 L ft 6 w ft 4 r WJ a(9) 

Let us now split the above integral in two parts: 

fill 8m 3 f 27T 
h = ~m (2^31^/0 ^«(*) 0* <*»(*) ^M*) -«<*»(*)) 

2m 2 f 27T 

d9(a 2 + B 2 -4cos{9)sm(9)a8)(8cos(9)sm(9)- a cos 2 (9)) 



37T 2 ft 5 Ju 

|"Ja( a 2 + 2/? 2 ), (110) 

h 1 1 19m 2 f 2rr 

'""' £ F / d6»(/3cos(6»)sm(0)-acos 2 (6»)) 



m (2tt) 2 3 ft 4 * 7 
— ^a. (Ill) 
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The spin-current T X y is thus given by 

% v « — ^a a 2 + 2/3 2 + — -f a. 112 

Note that the term cubic in SOI agrees with earlier results [5, 7], while the linear term survives here due to the absence 
of an anomalous velocity term. However, since ma(9) 2 <C h 2 Ep the cubic term is negligible and we obtain 



mEp 
irh 3 



% y » -^-a- (113) 



A similar calculation shows that 



r w «-^. (ii4) 

Note that differently from the main text, both Hamiltonian (100) and the spin current density are expressed in terms 
of Pauli matrices and not in terms of spin operators. Therefore, we need to multiply our results by a factor (H/2) 2 . 
Another multiplication with a factor h arises from the fact that the spin-orbit part is multiplied by 1/H in (100) as 
compared to our definition in the main text, 



mEp 

47T 



Txy « — — "> (115) 



r vy « -^0. die) 

VIII. SPIN-ORBIT INTERACTION WITH a = /3 

Let us consider the special case when Rashba and Dresselhaus coefficients are equal, i.e., a = ft. In such a case, 
the spin-orbit Hamiltonian takes the following form 

Hso = aY t {p B i +ft){Si-n\ ( 117 ) 

i 

and we define a gauge transformation U = e*^fc A *' r * with gauge vector field A& = (— am(s k — s v k ), —am(s k — s v k ), 0). 
Since 

Up^U' 1 = e iE * A *-**^ ,B e- <E * A *-**=^' , ' + om(«f-*V). ( 118 ) 

we have that 

Lr y Pi Lr -i = y (Pf + ""*(«? - a")) 2 + (g? + am ( j f - tf)) 2 = y gf , g , " 2 ™ / ng x 

•^ 2m ^—^ 2m / -—' 2m 2 



where the constant a 2 m/2 can be neglected without loss of generality. Under gauge transformation the spin isotropic 

•'3- 



term Hj = J^/j Sj • Ij transforms as follows, 



N T JVj JVj AT e 

UHjlf- 1 = UJ^Sj-IjU- 1 = jJ^USjU- 1 -Ij = jY j e i ^ Ak - H Y j s i 6{r i -R j )e- i ^ Ak - ik i j 

j=l j=l j=l »=1 

JVj JV e AT/ iV e 

= JEE . eiEfcAfc " fc§ f" tEfcAfc ' ffc , ' 5 ^ - RjA = J E E s ^ - R. WRi)Si • ij 

7=1 i— 1 „,„ x. 7 — 1 i=l 

AT 7 JV e Nr N e 

= jE E *(*< - R j) § ' • ^(Rj) Tf J = J T,J2 s & - R j) § ' • if 7 "- ( 12 °) 
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where 1Z <G 5*0(3) is a 3 x 3 special orthogonal matrix. The lattice spin Hamiltonian H% remains unchanged under 
gauge transformation. From this, we can conclude that Hamiltonian H = H e + H% + Hj + i?so & n d Hamiltonian 

i7KT = E^+E t/ ^-^) + E t/ ^)+E^- i '+ J E^-i? T J (i2i) 

i i<j i j,l j=l 



=H e =H X 

are equivalent, i.e., C/ _1 _ff K U = H . Let us now introduce an impurity spin rotation U such that U~ 1 IjU — if : 

U- 1 H U = H nT , (122) 

where Hq = H e + Hi + X^,=i &j ' Ij- We conclude that 

U- l U- l H UU = H. (123) 

The canonical ensemble average of an operator O is given by {0)h = Tr(e~* H 0)/Tr(e~P H ). Therefore, if we use the 
fact that U^U^HqUU = H, we obtain 

(0) H = n (e-&UHU-iu-* \uuOU~W~ 1 ]) /Tr(e~ fiUuHU - 1U - 1 ) =Tr(e- f>H <> [uUOU^lJ- 1 }) /Tr ( e -P H «) 
= (UUOU-W- 1 )^. (124) 

For the impurity spin magnetization fhi Eq. (124) implies that 

{rra) H = Tr (e" w <) /Tr (e"^) - (m?>* , (125) 

where to^ = UUrhiU^ 1 !/^ 1 = UrfiiU^ 1 is defined analogously as rhj but with rotated spins If = TZ(Rj)Ij. The 
explicit form of the matrix 1Z(Rj) is given by 

/ l/2 + l/2cos(0j) -l/2 + l/2cos(0j) -l/\/2sin(^)\ 
ft(Rj) = -1/2 + 1/2 cos((9j) l/2 + l/2cos(6> j ) -l/v/2sin(0j-) , (126) 

V l/\/2sin(0j) l/\/2sin(^) cos(^) / 

and it corresponds to a rotation around the axis n = (l/\/2, — 1/a/2, 0) by an angle of 9j = —\/2am(Rj + R v -). Note 
that the lattice vectors Rj = (R? , RX , Rf) should not be confused with the rotation matrix 1Z. 



A. (Anti-) ferromagnetic ordering 

The aim of this section is to rule out ferromagnetic ordering in a system defined by Hamiltonian H = H e +Hz + Hj + 
Hso+Hz(Q), where H§q is defined as in Eq. (117) and the Zeeman term is chosen to be -Hz(Q) = h J2jli( e ~ ' Rj Ij + 
e iQ-Rjjzy The magnetization is given by TO |(Q) = (mf(Q)) H = ^j(J2fli( e ' lQ ' Rj Ij + e iQ - R jif)) H - We know 
that fff(Q) = UUH z {Q)U- l lJ- 1 = /iEj2i( e_iQ ' Rj '(-^ f + e lQ '' R ' 3 {if) z ) and from Eq. (125) that (mj(Q)) H = 
((mfy(Q)) Ho = ^(E?±i(e- iQ ' Ri (ifY + e iQ ' Rj {lfY))H , with H a = H e + H T + Hj + Hf{Q). Since the rotation 
matrix 72. corresponds to a rotation around axis n = (l/\/2, — 1/\/2, 0) (see Eq. (126)), we define a new coordinate 
system, namely: z -> z' = (l/>/2, -l/\/2,0), y -)• y' = (l/\/2, l/\/2,0), and x ->• x' = z = (0,0,1). In this new 
coordinate system we obtain 

fff (Q) = - V (V i(Q ' +Q) ' Rj /+' + e i(Q'-Q)-Rj /-' + e i(Q-Q')-Rj /+' + e »(Q'+Q)-R i / 7 ' , j j ( 12 7) 

((m? ) z (Q))h - — ( E (e-W +Q >M/ + e i(Q'-Q)-Ri/-' + e *(Q-Q')-R 3 /+' + e i(Q'+Q)-Rj /-') \ ,(128) 

2Nl W ' /if 



where I* — If ± il y , and Q' = (— \/2am, — \/2am, 0). Note that in the new coordinate system x'y'z', Hz and i? 



J 



does not change their form. 
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The proof in section IV can be generalized in a straightforward way to a slightly different choice of operator 
Aq = /q +K — -t k + ^q+K' — ^q-K' an d symmetry-breaking Zemman term iJ^(K,K') = hJ2jli( e ~ tK Rj if + 
giK-R,- f-' + e -iK'.Rj /+' + giK'-Rj/-'^ with K = Q + Q' and K' = Q' - Q. We thus conclude that in the limit of 

vanishing external magnetic field (h — > 0) ((rhf) z ) u a = and consequently that (tHj)h — 0. Since the proof is valid 
for any possible Q, we have ruled out (anti-) ferromagnetic ordering along the z-direction in this case, too. 

B. Helical Ordering 

Here we rigorously rule out helical ordering in the direction of n = (1/V2, — l/v2, 0) for a system described 
by the Hamiltonian H = H e + H% + Hj + Hso + Hl£ , where Hgo is defined as in Eq. (117) and the Zeeman 
term is H^ = h^2jl 1 (If' ) z ■ The magnetization is given by (mf ) z = {{rhf ) z )h = ]\r(Sj=i(-^ Y)h- From 
UUHTflJ- 1 !)- 1 =H Z = hY^ii Ij and Eq. (125), we have 

((mf T y)H = {m}) Ho , (129) 

where Hq = H e + Hz + Hj + Hz- From the proof in section III, we know that in in the limit of vanishing external 
magnetic field (h —> 0) (itiJ)h = 0. From Eq. (129) we can thus conclude that ((mf' ) z )h = for h —J- 0. In order 
to have a better idea of the type of helical order excluded here, we can rewrite the magnetization (mf' ) z in the 
transformed coordinate system x'y'z' defined in the paragraph above, 

(<> = ^/f;(e-^M;'w*M-')\ , (130) 

1 " = 1 ' H 

where Q = (\/2am, \/2am, 0). This corresponds to a helix in the rr'y'-planc. 
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